Polariton Mott insulator with trapped ions or circuit QED 
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We consider variants of the Jaynes-Cummings-Hubbard model of lattice polaritons, taking into 
account next-nearest-neighbor, diagonal and long-range photon hopping in one and two dimensions. 
These models are relevant for potential experimental realizations of polariton Mott insulators based 
on trapped ions or microwave stripline resonators. We obtain the Mott-superfluid phase bound- 
ary and calculate excitation spectra in the Mott phase using numerical and analytical methods. 
Including the additional hopping terms leads to a larger Mott phase in the case of trapped ions, 
and to a smaller Mott phase in the case of stripline resonators, compared to the original model 
with nearest-neighbor hopping only. The critical hopping for the transition changes by up to about 
50 percent in one dimension, and by up to about 20 percent in two dimensions. In contrast, the 
excitation spectra remain largely unaffected. 

PACS numbers: 71.36. -fc, 73.43.Nq, 78.20.Bh, 42.50.Ct 



I. INTRODUCTION 

The study of strongly correlated and condensed phases 
of ultra-cold atoms has become an extremely active and 
versatile research field over the last decade. It has been 
triggered by the experimental realization of a Mott insu- 
lator to superfluid quantum phase transition of neutral 
atoms in an optical lattice ^ as described by the seminal 
Bose- Hubbard model [5] . The recent realization of strong 
hght-matter interaction in various cavity-quantum elec- 
trodynamics (QED) systems has triggered an immense 
interest in realizing strongly correlated and condensed 
phases with photons as well. Bose-Einstein condensation 
and superfluidity of weakly interacting polaritons, i.e., 
quasiparticles that form when photons strongly interact 
with matter, have already been observed experimentally 
with exciton-polaritons in a quantum well coupled to a 
photonic crystal cavity |3]. 

Today, a key challenge is to reach the limit of strong 
correlations, that is a photonic or polaritonic Mott in- 
sulator [IHS]. The Jaynes-Cummings-Hubbard model 
(JCHM) has been introduced to describe such an ex- 
treme state of light in an array of coupled high-Q elec- 
tromagnetic resonators In the generic JCHM, pho- 
tons hop between nearest-neighbor (nn) cavities and in- 
teract with a single two-level system locally in each cav- 
ity. This fundamental light-matter coupling introduces 
a non-linearity into the system, leading to an effective 
repulsive photon-photon interaction. The competition 
between hopping (delocalization) and light-matter inter- 
action (localization) leads to a phase diagram featuring 
Mott lobes similar to the Bose-Hubbard model. The 
phase diagram, excitations and critical exponents of the 
generic JCHM have been studied theoretically using var- 
ious numerical and analytical methods fiWTF and it has 
been shown that the same critical theory applies to the 
JCHM and Bose-Hubbard model [TOHni E]- Original 
proposals for an experimental realization of the JCHM 



were based on nitrogen vacancy centers in diamond [1] 
or self-assembled quantum dots in photonic crystals |16| . 
However, both variants suffer from short coherence times 
and are prone to disorder effects. 

Two rather recent proposals, which are not limited 
with respect to decoherence and disorder, are based on 
circuit QED [TUIIIZIIIH] and trapped ions [Ullin]. In cir- 
cuit QED, Josephson qubits are coupled to stripline res- 
onators. Arrays can be formed by capacitively coupling 
resonators in one-dimensional (ID) or two-dimensional 
(2D) geometries. In trapped ions, the role of photons is 
played by phonon excitations which are exchanged be- 
tween the ions, and arrays can be formed either in linear 
Paul traps or in microtraps. However, in both cases, pho- 
ton transfer beyond the canonical nn hopping cannot be 
neglected. In a circuit QED setup, the 2D square lat- 
tice geometry gives rise to next-nearest neighbor (nnn) 
hopping. In the case of trapped ions, hopping is medi- 
ated by the Coulomb interaction between dipoles, and 
is thus long ranged in all lattice dimensions. An accu- 
rate quantitative description of the Mott-superfluid tran- 
sition in these systems requires more general versions 
of the JCHM. In this paper, we show that additional 
hopping terms substantially modify the Mott-superfluid 
phase boundary. We compare exact numerical results 
obtained from extensive quantum Monte Carlo (QMC) 
simulations with the variational cluster approach (VCA) 
and analytical calculations based on a linked-cluster ex- 
pansion. 

The paper is organized as follows. In Sec. |TI| we intro- 
duce the JHCM with a general hopping term, and specify 
the cases pertaining to circuit QED arrays and trapped 
ions. In Sec. |IIH we outline the numerical and analytical 
methods used. Section IIVI contains a discussion of our 
results for the effect of additional hopping terms on the 
phase diagram and the elementary excitations. A sum- 
mary is given in Sec. |V] Finally, the Appendix contains 
relevant results for the atomic limit. 
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FIG. 1: (Color online) Circuit-QED realization of the JCHM 
on a 2D square lattice. A filled circle indicates a lattice site, 
defined by the position of a qubit. Each qubit is located inside 
a resonator represented by a rectangle. The capacitive cou- 
pling between resonators gives rise to photon hopping across 
the lattice. Two different couplings occur: nearest-neighbor 
hopping with amplitude t along the lattice axes, and diagonal 
next-nearest neighbor hopping in every other plaquette with 
amplitude t' . 



II. MODELS 

The JCHM is defined by the Hamiltonian 

^ = E + E ^^M'^J + H.c.) -fiN, (1) 

i ij 

with the on-site JC Hamiltonian for the cavity at site i 
of the lattice, 

ftJ^ = Wboja^ + UJ^atcr' + 5(^+0, + a^^ aj) . (2) 



A. Trapped ions 

The case of trapped ions implies a ID model with frus- 
trated long-range hopping 



arising from dipole-dipole interactions between confined 
ions [T9j. Here, t denotes the nn hopping strength. The 
effects of long-range hopping for trapped ions have been 
investigated using simple mean-field Hke approximations 
|20) . However, the validity of mean- field theory in such a 
ID setting is not clear. In particular, it fails to predict the 
Kosterlitz-Thouless transition at the tip of the Mott lobe 
which manifests itself as a non-analytic cusp in the phase 
boundary [71[2I]. In this paper, we therefore address this 
problem using numerical methods. 

Due to the alternating sign in Eq. ([3|, which gives rise 
to frustration, QMC simulations are impracticable. In- 
stead, we use the VCA, a quantum cluster method which 
can be applied to frustrated systems and also to com- 
pute single-particle excitation spectra. To benchmark the 
VCA for the case of long-range hopping, we also present 
results for the case 

i.e., without the alternating sign. Hamiltonian ([T]) with 
the hopping integrals given by Eq. Q can be studied us- 
ing QMC simulations to obtain the exact phase boundary 
for the Mott-superfluid transition. 



Here, Wb denotes the energy of bosons (photons or 
phonons) , and Wq is the level splitting of the qubits (ions 
or Josephson qubits). The spin operators ct* describe in- 
trasite transitions between the two qubit levels induced 
by emission or absorption of a boson with rate g. The 
bosonic operators a, fulfill the usual commutation re- 
lations [tti, Oj] = Sij. The coupling g gives rise to the for- 
mation of polaritons (combined boson-qubit excitations) 
whose number operator N = '^i{aiai+(T^iJ~) commutes 
with H. Hence, within this theoretical framework, the 
polariton number is conserved and can be controlled via 
the chemical potential fi. 

The second term in Eq. ([T]) describes photon hopping 
between site i and site j of the cavity array with ampli- 
tude tij. At this point, the hopping is completely gen- 
eral, and may involve any pair of lattice sites. In previous 
work, only nn transfer along the directions of the lattice 
basis vectors was considered, with — —t. We refer to 
this case as the original JCHM in the following. 

Here we consider the JCHM with more complex hop- 
ping terms, motivated by possible realizations using 
trapped ions in ID Paul traps, and stripline resonators 
on a 2D square lattice. 



B. Stripline resonators 

The 2D JCHM may be realized in experiment using 
stripline resonators on a square lattice. A simple config- 
uration is shown in Fig.[T| The hopping integral = —t 
when i and j are nn's along a lattice bond (as in previ- 
ous studies of the 2D JCHM), and tij = —t' for diagonal 
nnn's on every other plaquette. For any other pair of 
lattice sites we have tij = 0. 

In this circuit QED setup, the hopping rates t and t' 
are proportional to the mutual capacitances of the res- 
onators in Fig. [T] [171 ES] • By choosing a proper geometry 
and/or distance between pairs of resonators on the diag- 
onal of the array (determining t') and pairs of resonators 
on the lattice axes (determining t) , the diagonal hopping 
rate t' can be made much smaller with respect to the nn 
hopping rate t and vice versa. Thus, the ratio f /t can be 
engineered almost arbitrarily [23] . Here we consider the 
ratio i'/i = 1/2. 

The freedom to tune the ratio t' /t to substantial values 
in circuit QED systems is in strong contrast to ultra-cold 
atoms in optical lattices. For the latter, hopping inte- 
grals beyond nn pairs are orders of magnitude smaller 
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[21]. Moreover, diagonal hopping is absent on hypercu- 
bic lattices because of the orthogonality of the Wannier 
states. The small size of hopping integrals beyond the nn 
terms makes their effect very small, with correspondingly 
few studies available . 

The circuit QED model with hopping integrals t and 
t' as depicted in Fig. [l] can be simulated using the QMC 
method. For comparison, we also present results for a 2D 
model with nn hopping t, and nnn hopping t" between 
sites separated by two lattice constants, along the (1,0) 
and (0, 1) directions. 



III. METHODS 

This section describes the methods used in this work 
in some detail. In two dimensions, we use exact QMC 
simulations, the VGA, and an analytical linked-cluster 
approximation around the mean-field limit in two dimen- 
sions. Since mean-field theory is not a valid starting point 
for an expansion in one dimension, we show only numer- 
ical results (QMC, VCA) for that case. 



A. Quantum Monte Carlo 

In the absence of frustration, bosonic Hamiltonians of 
the form of ([I]) can be studied by exact QMC simulations. 
A particularly popular representation is the continuous- 
time stochastic series expansion [2B ; it has been applied 
to the JCHM model before i9iLl5, 27J. We use the ALPS 
1.3 implementation [28] of the stochastic series expansion 
with directed loop updates [2HHSI]- A maximum site oc- 
cupation of six photons is sufficient to make the trun- 
cation error negligible. The long-range hopping defined 
by Eq. Q is treated by taking into account all transfer 
processes with \i — j\ < L/2 — 1 for system size L (due to 
periodic boundary conditions). 

The Mott-superfluid phase boundary can be deter- 
mined by calculating the superfluid density 

{^') 

where D denotes the dimension of the lattice, w is the 
winding number |32j . and /3 is the inverse temperature. 
Exploiting the scaling form of ps j2|, the critical point 
can be obtained from simulations with different system 
sizes at inverse temperature P/L^ = const. The dynami- 
cal critical exponent of the JCHM is z = 1 for the fixed- 
density transition, and z = 2 for the generic transition 
[15j . For large enough values of the system size L, curves 
for different L values intersect at a single point [HI [H] 
which defines tc (for horizontal scans in the fi phase dia- 
gram) respectively (for vertical scans) . We use system 
sizes up to L = 64 in one dimension, and up to 40 x 40 for 
the 2D models. The inverse temperatures were typically 
/3/L2 = 1/4 for z = 2, and P/L = 2 or 4 for z = 1. The 



resulting accuracy for the phase boundary in units of g 
is estimated to be better than 0.0005 (below the symbol 
size used in the figures). 

B. Variational cluster approach 

The minus-sign problem resulting from the frustration 
induced by choice Q for the hopping integrals motivates 
the use of an alternative numerical method. The VCA 
is a cluster method; hopping within a finite reference 
cluster is treated exactly, and hopping beyond the clus- 
ter is taken into account perturbatively. For the relation 
to other popular cluster methods such as the dynamical 
mean field theory, see [3^. The VCA [53] can be ap- 
plied for any choice of hopping integrals in one and two 
dimensions and, also, permits calculation of excitation 
spectra. It has been applied to the original JCHM with 

= in 8J and |33j , where in the latter case the problem 
was mapped to an effective model. In one dimension, the 
VCA yields quantitatively reliable results f51fM|. In 2D, 
accurate results can be obtained except for a small region 
around the Mott lobe tip in the phase diagram. 

The quality of the approximation depends on the size 
of the cluster compared to the correlation lengths of 
the problem, and the number of variational parameters. 
Since larger clusters and more parameters have compa- 
rable impact, we use only the boson energy in the cluster 
as a parameter, and vary the cluster size. This is also 
motivated by the presence of long-range hopping terms 
in the models considered here. For more details see [5]. 

Here we investigate Hamiltonian ([!]) using the VCA in 
its formulation for bosonic systems 35,. We work exclu- 
sively at zero temperature. To deal with the long-range 
hopping defined by Eqs. ([s]) and Q, we allow for hopping 
processes up to a distance L — 1, where L is the cluster 
size used in the calculation. 

On the technical side, we note that care is required 
when carrying out the sums over wave vectors in the cal- 
culation of the grand potential [8", '33] . Far from the lobe 
tip, the sums converge rapidly, whereas close to the tip 
an increasingly finer mesh is required. We attribute this 
effect to the momentum dependence of the self-energy 
and excitations. Deep in the Mott phase, the particle 
and hole bands are almost completely flat, whereas close 
to the lobe tip they acquire a substantial dispersion [5|. 
Given convergence, the VGA yields the correct form of 
the phase boundaries, leading to improved results for the 
2D case compared to exact QMC data [SI 12] . 

C. Linked-cluster expansion 

We analytically calculate the photonic Matsubara 
Green's function Gy(T;r') = —{Tai{T)aj{T')) with the 

time-ordering operator T and aj{T') = e^'^ aJjC^^'^ 
using a linked-cluster expansion in terms of local cu- 
mulants originally developed by Metzner [36j for the 
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Fermi-Hubbard model, and recently applied to the Bose- 
Hubbard model ^ and the JCHM ^12]. After a Fourier 
transformation and analytic continuation, the inverse 
Green's function directly yields the excitation spectrum 
a;(k) via G^^(k, = and the phase boundary tc{n) via 
G^^(0, O)lt^(^) = 0. Key results are summarized below. 



1. Random phase approximation 

Within the mean-field random phase approximation 
(RPA), the full Green's function is given by [TT] 



with the local one-particle cumulant 



(6) 



(7) 



The hole/particle spectral weights and bare excitation 
energies in the numerator and denominator in Eq. ([t]) 
are given in the Appendix. In Eq. ([6|, J(k;) denotes the 
bare dispersion. For the 2D circuit QED model in Fig. [l] 



J(k) — 2t(cos kx + cos ky 



2t' cos kxky . 



(8) 



From Eq. ([6]), we can derive analytical expressions for 
the phase boundary and the excitation spectrum. For 
example, the RPA phase boundary for the t, t' model at 
a fixed ratio R = t'/t is given by l/U = 4G°(0)(1 + 
R/2). Thus, on the mean-field level, nnn hopping simply 
rescales the nn hopping according to tc ic/(l + R/2). 
Hence, a positive (negative) nnn hopping decreases (in- 
creases) the size of the Mott lobes and makes mean- 
field theory a better (worse) approximation. The RPA 
is equivalent to the cluster perturbation theory (VGA 
without variational parameters) for a single-site cluster 
[35|. 



2. One-loop approximation 

The mean-field RPA corresponds to a summation of 
all self-avoiding walks through the lattice. In order to 
take into account the leading quantum correction one 
can additionally include all one-time forward/backward 
hopping processes between two neighboring sites. The 
resulting so-called one-loop approximation was previously 
found to be in excellent agreement with numerical results 
for the phase diagram for both the Bose-Hubbard model 
[37j and the original JCHM with = [HJ- 

IV. RESULTS 

To calculate the phase diagram of the models discussed 
in Sec. [ill we vary the nn hopping strength t common to 



all models, and scale any additional hopping integrals 
accordingly. In the following, we use g as the unit of 
energy, and consider Wb = = 1 (resonance condition). 
We also set h, /cb and the lattice constant to 1. 



A. Mott-superfluid transition 

Similarly to cold atoms in optical lattices, the JCHM 
can be tuned across the Mott-superfluid transition by 
changing the ratio t/g. We therefore expect the addi- 
tional hopping terms discussed in Sec. |ll] to modify the 
extent of the Mott insulating region in the phase dia- 
gram. Because the Mott phase is of particular interest 
as an initial state for possible applications in quantum 
computing and quantum information, we determine the 
size of the first (largest) Mott lobe with polariton density 
n — I. 



1. One dimension 

We begin with the ID model with long-range hopping 
given by Eq. Q, which represents an interesting theo- 
retical problem for two reasons. First, it allows us to 
benchmark the VGA against exact QMC results for the 
novel case of long-range hopping. Second, the long-range 
hopping is expected to mimic the effect of increasing the 
lattice coordination Z , while keeping the lattice topol- 
ogy and dimension unchanged. Previous VGA calcula- 
tions for the JGHM have demonstrated a crossover from 
Kosterlitz-Thouless behavior (refiected in a strongly non- 
parabolic shape of the phase boundary, including reen- 
trant behavior [^) on a ID chain, to a mean-field like, 
parabolic phase boundary for the 2D square lattice [5]. 
This crossover is also visible upon comparing Figs. [2] 
and|3] Mean-field behavior in a ID matter-light system 
is also suggested by results for a two-component Bose- 
Hubbard model coupled to a global photon mode (as in 
the Dicke model) [38j; the coupling gives rise to photon- 
mediated long-range interaction. 

Figure [2|^a) shows the Mott-superfluid phase bound- 
ary for the first Mott lobe. The long-range hopping sub- 
stantially reduces the extent of the insulating phase, as 
compared to the original JGHM (shaded regions, results 
taken from T ). Due to the particular shape of the phase 
boundary in one dimension, the reduction of the critical 
hopping due to tij can be very large, about 50 percent 
for fji/g — 0.1. 

The Kosterlitz-Thouless phase transition at a fixed po- 
lariton density, which occurs at the tip of the Mott lobe, 
is difficult to study by the QMG method. Away from the 
lobe tip, the critical points from QMG agree well with 
the VGA results. The approximate cluster approach un- 
derestimates the effect of quantum fluctuations for larger 
i/g, leading to a slightly larger Mott phase than found 
by QMG as visible in Fig. ^a) for the QMG data points 
located at t/g > 0.1. This flnding is in accordance with 
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FIG. 2: (Color online) Zero-temperature Mott lobe with den- 
sity n = 1 of the ID JCHM, showing the effect of (a) long- 
range hopping defined by Eq. (j4|, and (b) sign-alternating 
long-range hopping defined by Eq. Q, as appropriate for an 
experimental realization based on trapped ions. Case (a) can 
be investigated by QMC simulations (symbols), and the phase 
boundary is compared to the VCA (lines, for a cluster size 
L = 8). For reference, we also show the exact result for 
the original JCHM with nearest-neighbor hopping only [7] 
(shaded regions). Case (b) is not accessible for QMC, and we 
rely on the VCA for the phase diagram. 
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FIG. 3: (Color online) Zero-temperature Mott lobe with den- 
sity n = 1 of the 2D JCHM, showing the effect of (a) next- 
nearest-neighbor hopping t" — t/8 (see Sect. [ll| and (b) di- 
agonal next-nearest-neighbor hopping t' = t/2 as appropriate 
for an experimental realization based on stripline resonators 
(see Fig. [l] and Sec. [ll|. We show results from QMC simu- 
lations (symbols), the VCA (lines, L = 8), the mean-field 
theory (RPA), and the one-loop approximation. For refer- 
ence, we also include the t' = t' = Q QMC results for the 
original 2D JCHM model (shaded regions, taken from ^). 



the cluster size dependence of the VCA Mott lobe in 8J, 
where satisfactory convergence was achieved for L = 8 
(the same size as used here). Far away from the lobe tip, 
for t/g < 0.1 in Fig. [2|a), the transition is essentially 
driven by density fluctuations, and the VCA and QMC 
results agree almost perfectly. 

Concerning a possible change of the universality class 
in the presence of long-range hopping, neither the QMC 
nor the VCA results provide any evidence of such a sce- 
nario. The shape of the Mott lobe for the model with 
long-range hopping is reminiscent of the phase bound- 
ary of the original JCHM (shaded region). Besides, the 
QMC results for the superfluid density (not shown) ex- 
hibit scaling with dynamical critical exponent z = 1. The 
absence of mean-field behavior despite long-range hop- 
ping is a consequence of the fast decay with distance, 
tij ~ |»- j^^ see Eq. ^. 



Turning to the model for trapped ions, with sign- 
alternating long-range hopping as given by Eq. ([s]), we 
note that the QMC method suffers from a severe minus- 
sign problem. Relying on the accuracy of the VCA es- 
tablished in Fig. ^a) and in previous work, we show the 
VCA phase boundary in Fig. . The sign-alternating 
hopping has an effect opposite to the choice ([4|. The 
extent of the Mott is noticeably increased compared to 
the model with nn hopping only. The largest change of 
the critical hopping again occurs for fi/g w 0.1, with an 
increase of about 50 percent. The qualitative shape of 
the phase boundary is unchanged. 



2. Two dimensions 

We now turn to the results for the 2D models discussed 
in Secjll] Since there is no frustration in either case, we 
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FIG. 4: (Color online) Zero-temperature excitation spectra 
of the ID JCHM at zero detuning with (a) nearest-neighbor 
hopping t, and (b) sign-alternating long-range hopping Uj, 
relevant for trapped ions [Eq. (|3|] . Results are from the VGA 
with cluster size L = 8, and using t/ g = 0.1, ^/g — 0.15, see 
Fig.^b). 



can study the Mott-superfluid transition exactly using 
the QMC method. The findings are compared to the 
VGA, and to analytical approximations. 

Figure [sj^a) contains the results for the t, t" model, 
with t" /t =1/8 [motivated by taking a.\/r^ dependence 
as in Eq. Q but keeping only the first two terms]. The 
hopping t causes a 10 percent reduction of the critical 
hopping close to the tip as compared to the model with 
nn hopping only [9]. The VGA provides a remarkably 
good description of the whole phase boundary. In partic- 
ular, the spurious inflection points visible in earlier work 
[HI |5] are absent. As in one dimension, the underestima- 
tion of spatial fluctuations leads to slightly larger values 
of the critical hopping for the Mott-superfluid transition. 
Gonversely, the RPA mean-field result severely underes- 
timates the effect of fluctuations, a deficiency which is to 
a large extent remedied within the one-loop approxima- 
tion. 

Results for the 2D circuit QED model with hoppings 
t' /t = 1/2 are shown in Fig. [sjjb). The effect of t' is 
qualitatively very similar to that in Fig. [3]^a), with the 
larger value of t' compared to t" leading to a stronger 
decrease of the critical hopping close to the tip (about 20 
percent). The agreement between the different methods, 
in particular the VGA and QMC, is even better than for 
the t,t" model. Hence, except for details such as criti- 
cal exponents [15], both the VGA and the one-loop ap- 
proximation may be used for semi-quantitative analyses 
of the 2D JGHM at substantially smaller computational 
cost than for the QMG. 



FIG. 5: (Golor online) Zero-temperature excitation spectra 
for the 2D JGHM at zero detuning with t' = Q (solid lines), 
or t' /t = 1/2 [dashed lines and intensity plot in (a)] as appro- 
priate for a circuit-QED setup (see Fig.[T]and Sec.|ll|. (a) Pa- 
rameters deep in the Mott lobe, using t/g = 0.02, /I/p = 0.22. 
(b) Parameters exactly at the Mott lobe tip, using t/g = 0.04 
(for t' = 0) respectively t/g = 0.032 (for t' /t = 1/2), and 
^,/g = 0.22, see Fig. |3][b). Lines are RPA results for the 
particle/hole dispersion. The density plot in (a) shows VGA 
results based on a 2 x 2 cluster. 



B. Excitation spectra 

Single-particle excitation spectra are directly accessi- 
ble in cavity QED via photocmission spectroscopy. As 
demonstrated previously ^JJ, 27j, they allow us to iden- 
tify the system state via the presence or absence of a Mott 
gap and to calculate effective particle and hole masses. 
We therefore study here the effect of the additional hop- 
ping processes in models for trapped-ion and circuit-QED 
realizations of the JGHM. 

Figure |4] compares the excitation spectrum deep in 
the Mott phase for the cases of the original ID JGHM, 
Fig. |4][a), and the model with sign-alternating long- 
range hopping, Fig. |4][b). Results were obtained us- 
ing the VGA. The four gapped branches of the original 
model persist in the presence of long-range hopping; the 
two low-lying branches correspond to conventional parti- 
cle/hole excitations [5], whereas the high-energy branches 
are so-called conversion modes, which arise from the com- 
posite nature of polaritons [11] . The upper modes have a 
very small bandwidth and low spectral weight. Gompar- 
ing Figs. |4][a) and (b), we conclude that the effect of the 
long-range hopping is negligible. A comparison close to 
the lobe tip is difficult in one dimension, because of the 
different extent of the Mott phases, see Fig.|2j[b). 

Results for the excitation spectra of the 2D circuit- 
QED model are shown in Fig. [s] focusing on low-energy 
modes. Deep in the Mott phase, we find perfect agree- 
ment between the VGA and the RPA. The renormaliza- 
tion of the Mott-superfluid transition (including the Mott 
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gap) due to t' [see Fig. [3|Jb)] is visible from the particle 
and hole bands close to the T point. 

Finally, we consider the spectrum exactly at the lobe 
tip. Because our formulation of the VGA is restricted to 
the gapped Mott phase, we only show the corresponding 
RPA results in Fig. Sjb). The ratio t/g has been tuned 
to the lobe tip for both i' = and t' > 0. The RPA 
dispersion becomes relativistic with gapless, linear modes 
in both cases, and noticeable differences due to t' only 
far away form the F point and hence at high energies. In 
particular, the sound velocities are practically identical. 
The presence of linear excitations at the lobe tip is a 
signature of the fixed-density transition with dynamical 
critical exponent z = 1. The value z = 1 has recently 
been demonstrated using QMC simulations [151 IM] for 
the case t" = t' = (original JCHM). As expected, the 
universality does not change in the presence of additional 
hopping terms. 



V. SUMMARY 

We have investigated the effects of additional hopping 
terms beyond the generic nearest-neighbor transfer in the 
Jaynes-Cummings-Hubbard model in one and two dimen- 
sions. Such hopping terms arise when considering possi- 
ble experimental realizations based on trapped ions in a 
linear Paul trap, or stripline resonators on a square lat- 
tice. Using numerical and analytical methods, we have 
shown that the phase diagram is modified substantially. 
Compared to the original model with nearest-neighbor 
hopping only, the Mott lobe with density 1 becomes en- 
larged in the case of trapped ions, but is reduced for 
the case of stripline resonators. This effect is particu- 
larly pronounced at selected chemical potentials in one 
dimension, due to the shape of the Mott lobe. In con- 
trast, excitation spectra are only very weakly affected by 
the additional hopping terms, especially at low energies. 
In particular, the sound velocity at the lobe tip as well 
as the universality class of the JCHM remain unchanged. 

An interesting problem arises from the possibility of 
tuning the hopping ratio t' /t in a circuit-QED setup. For 
t' ^ t, the lattice separates into independent, diagonal, 
one-dimensional chains. By varying t' /t, one could thus 
observe a crossover from one to two dimensions, i.e., from 
Kosterlitz-Thouless behavior to mean-field like behavior. 
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Appendix: Atomic-limit results 

In the atomic limit {Uj = 0) the cigcnstatcs of the 
Hamiltonian ([T]) are the dressed polariton states labeled 
by the polariton number n and upper/lower branch index 
CT = ±. For n > they can be written as a superposition 
of a Fock state with n photons plus atomic ground state 
|n,g) and {n — 1) photons with the atom in its excited 
state \{n — 1), e), 

\n+) = sin6'„|?i,g) + cos^nKn - 1) ,e) , 

|n-) = cos6'„|n,g) - sin^nKn - 1) ,e) , (A.l) 



where 



with 



tan0„ = 2gy/n/{6 + 2x„) 



(A.2) 
(A.3) 



and the detuning 6 ~ uj]^ — coq. The eigenvalues are 

-i^i-iOh)n-d/2 + aXn, a ^ ± . (A.4) 



The zero polariton state |0- 



\0,g) is a special case 



with Cq ~ 0. Upper and lower polariton energies are 
separated by the Rabi splitting r2„ = 2x„. The spectral 
weights in Eq. ([T]) are then defined as 



(A.5) 



with the matrix elements f^''^ — {n(j\a^\{n — l)<j')- The 
bare excitation energies are given by 



(A.6) 
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